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We describe the antenna subtraction method for treating real emission singularities in the calculation of jet 
observables at NNLO accuracy, in particular in view of the computation e + e~ — > 3 jets at NNLO. 



1. Introduction 

Using experimental data on jet production ob- 
servables to extract QCD parameters, especially 
the QCD coupling constant a 8 , one finds that the 
dominant source of error on these parameters is 
very often the uncertainty inherent to the QCD 
calculations of the obervables. These calculations 
are at present available to next-to-leading (NLO) 
accuracy, their precision can be improved only by 
calculating corrections to the next perturbative 
order: next-to-next-to-leading order (NNLO). 

In the recent past, many ingredients to NNLO 
calculations of collider observables have been de- 
rived (see and references therein): the mass- 
less two-loop 2^2 and 1 — > 3 matrix elements 
relevant to NNLO jet production have been com- 
puted and are now available for many processes 
of phcnomenological relevance. The one-loop cor- 
rections to 2 — ► 3 and 1 — ► 4 matrix elements have 
been known for longer and form part of NLO cal- 
culations of the respective multi-jet observables. 
These NLO matrix elements naturally contribute 
to NNLO jet observables of lower multiplicity if 
one of the partons involved becomes unresolved 
(soft or collinear). In these cases, the infrared 
singular parts of the matrix elements need to be 
extracted and integrated over the phase space ap- 
propriate to the unresolved configuration to make 
the infrared pole structure explicit. As a final in- 
gredient, the tree level 2^4 and 1 -t 5 pro- 
cesses also contribute to (2 — » 2)- and (1 — > 3)- 
type jet observables at NNLO. These contain dou- 
ble real radiation singularities corresponding to 



two partons becoming simultaneously soft and/or 
collinear. To compute the contributions from sin- 
gle unresolved radiation at one-loop and double 
real radiation at tree level, one has to find sub- 
traction terms which coincide with the full matrix 
elements in the unresolved limits and are still suf- 
ficiently simple to be integrated analytically in 
order to cancel their infrared pole structure with 
the two-loop virtual contributions. It is the aim of 
this talk to present a systematic method, named 
antenna subtraction, to construct NNLO subtrac- 
tion terms and to illustrate its applications. 

2. Antenna subtraction 

An m-jet cross section at NLO is obtained by 
summing contributions from (m + l)-parton tree 
level and m-parton one-loop processes: 



da NLO = / {dcr NLO - do- 

J d$ m+ i 



5 ^ 
NLO) 



dfJ NLO 



dcr 7VLO 



The cross section da§ LO is the (m + l)-parton 
tree-level cross section, while da^ LO is the one- 
loop virtual correction to the m-parton Born cross 
section da B . Both contain infrared singularities, 
which are explicit poles in 1/e in da^ LO , while 
becoming explicit in dafj LO only after integration 
over the phase space. In general, this integration 
involves the (often iterative) definition of the jet 
observable, such that an analytic integration is 
not feasible (and also not appropriate). Instead, 
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Figure 1. Illustration of NLO antenna factori- 
sation representing the factorisation of both the 
squared matrix elements and the (to + l)-particle 
phase space. The term in square brackets repre- 
sents both the antenna function and the antenna 
phase space. 
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Figure 2. Illustration of NNLO antenna factori- 
sation representing the factorisation of both the 
squared matrix elements and the (to + 2)-particle 
phase space when the unresolved particles are 
colour connected. The term in square brackets 
represents antenna function and phase space. 



one would like to have a flexible method that can 
be easily adapted to different jet observables or 
jet definitions. Therefore, the infrared singulari- 
ties of the real radiation contributions should be 
extracted using infrared subtraction terms. One 
introduces da^ LO , which is a counter-term for 
d° nloi navm g the same unintegrated singular 
behaviour as da^ LO in all appropriate limits. 
Their difference is free of divergences and can be 
integrated over the (to + l)-parton phase space 
numerically. The subtraction term d<r^ LO has to 
be integrated analytically over all singular regions 
of the (to + l)-parton phase space. The result- 
ing cross section added to the virtual contribution 
yields an infrared finite result. Several methods 
for constructing NLO subtraction terms system- 



atically were proposed in the literature 
For some of these methods, extension to NNLO 
was discussed [Zj an d partly worked out. We focus 
on the antenna subtraction method |2I3| . which 
we extend to NNLO. 

The basic idea of the antenna subtraction ap- 
proach at NLO is to construct the subtraction 
term dcr^ LO from antenna functions. Each an- 
tenna function encapsulates all singular limits 
due to the emission of one unresolved parton be- 
tween two colour-connected hard partons (tree- 
level three-parton antenna function). This con- 
struction exploits the universal factorisation of 
phase space and squared matrix elements in all 
unresolved limits, depicted in Figure ^ The in- 
dividual antenna functions are obtained by nor- 
malising three-parton tree-level matrix elements 
to the corresponding two-parton tree-level matrix 
elements. 

At NNLO, the m-jet production is induced 
by final states containing up to (m + 2) par- 
tons, including the one-loop virtual corrections to 
(m + l)-parton final states. As at NLO, one has 
to introduce subtraction terms for the (to + 1)- 
and (to + 2)-parton contributions. Schematically 
the NNLO m-jet cross section reads, 



da 



NNLO 



d<I>„ 



(da 



R 

NNLO 



d*„ 



da 



NNLO ^ a NNLO ) 



d&NNLo) 



VS,1 \ 



d&NNLO 



da 



d*„ 



VS.l 

NNLO 



da 



V.2 

NNLO ' 



where da NNLO denotes the real radiation sub- 
traction term coinciding with the (to + 2)-parton 
tree level cross section dafj NLO in all singular 
limits 8 . Likewise, derJ^^Q is the one-loop vir- 
tual subtraction term coinciding with the one- 
loop (m + l)-parton cross section dcrJ^ io in all 
singular limits 9 . Finally, the two-loop correc- 
tion to the TO-parton cross section is denoted by 

aa NNLO- 

Both types of NNLO subtraction terms can be 
constructed from antenna functions. In dafj NLO , 
we have to distinguish four different types of un- 
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Figure 3. Illustration of NNLO antenna factorisation representing the factorisation of both the one- loop 
"squared" matrix elements (represented by the white blob) and the (m + l)-particle phase space when 
the unresolved particles are colour connected. 



resolved configurations: (a) One unresolved par- 
ton but the experimental observable selects only 
m jets; (b) Two colour-connected unresolved par- 
tons (colour-connected); (c) Two unresolved par- 
tons that are not colour connected but share a 
common radiator (almost colour-unconnected); 
(d) Two unresolved partons that are well sepa- 
rated from each other in the colour chain (colour- 
unconnected). Among those, configuration (a) is 
properly accounted for by a single tree-level three- 
parton antenna function like used already at 
NLO. Configuration (b) requires a tree- level four- 
parton antenna function (two unresolved partons 
emitted between a pair of hard partons) as shown 
in FigureEl while (c) and (d) are accounted for by 
products of two tree-level three-parton antenna 
functions. 

In single unresolved limits, the one-loop cross 



section da 



v,i 

NNLO 



is described by the sum of two 



terms 9 : a tree- level splitting function times a 
one-loop cross section and a one-loop splitting 
function times a tree-level cross section. Conse- 
quently, the one-loop single unresolved Subtrac- 
ts S 1 

tion term da N] ^ LO is constructed from tree-level 
and one-loop three-parton antenna functions, as 
sketched in Figure |2| Several other terms in 
^nnlo cance l with the results from the inte- 
gration of terms in the double real radiation sub- 
traction term dafj NLO over the phase space ap- 
propriate to one of the unresolved partons, thus 
ensuring the cancellation of all explicit infrared 



poles in the difference da 



v,l 

NNLO 



da 



Finally, all remaining terms in da 



vs.i 

NNLO- 



NNLO an< l 
V S 1 

da Nl ^ LO have to be integrated over the four- 
parton and three-parton antenna phase spaces. 



After integration, the infrared poles are rendered 
explicit and cancel with the infrared pole terms 
in the two-loop squared matrix element da 



V,2 
NNLO' 



3. Derivation of antenna functions 



The subtraction terms da NLO , da NNLO and 



da 



vs,i 

NNLO 



require three different types of antenna 



functions corresponding to the different pairs 
of hard partons forming the antenna: quark- 
antiquark, quark-gluon and gluon-gluon antenna 
functions. In the past |2l3j . NLO antenna func- 
tions were constructed by imposing definite prop- 
erties in all single unresolved limits (two collinear 
limits and one soft limit for each antenna). This 
procedure turns out to be impractical at NNLO, 
where each antenna function must have defi- 
nite behaviours in a large number of single and 
double unresolved limits. Instead, we derive 
these antenna functions in a systematic man- 
ner from physical matrix elements known to pos- 
sess the correct limits. The quark- antiquark an- 
tenna functions can be obtained directly from the 
e + e~ — * 2j real radiation corrections at NLO and 
NNLO ^U]. For quark-gluon and gluon-gluon an- 
tenna functions, effective Lagrangians are used to 
obtain tree-level processes yielding a quark-gluon 
or gluon-gluon final state. The antenna func- 
tions are then obtained from the real radiation 
corrections to these processes. Quark-gluon an- 
tenna functions were derived from the purely 
QCD (i.e. non-supersymmetric) NLO and NNLO 
corrections to the decay of a heavy neutralino 
into a massless gluino plus partons while 
gluon-gluon antenna functions |KSj result from the 
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QCD corrections to Higgs boson decay into par- 
tons [Tl) . 

All tree-level three-parton and four-parton an- 
tenna functions and three-parton one-loop an- 
tenna functions are listed in £Q, where we also 
integrate them, using the phase space integration 
techniques described in |15) . 

4. Application to e + e — > 3 jets 

To illustrate the application of antenna sub- 
traction on a non-trivial example, we derived 
in [TUB] the l/A 2 -contribution to the NNLO cor- 
rections to e + e~ — » 3 jets. This colour fac- 
tor receives contributions from 7* — > qqggg and 
7* ~^ at tree-level ^J], 7* ~~ * QQ99 an d 

7* — * QQQQ at one-loop [Hi] and 7* — > (795 at two- 
loops ^Hj- The four-parton and five-parton final 
states contain infrared singularities, which are ex- 
tracted using the antenna subtraction formalism. 

In this contribution, all gluons are effectively 
photon-like, and couple only to the quarks, but 
not to each other. Consequently, only quark- 
antiquark antenna functions appear in the con- 
struction of the subtraction terms. 

Starting from the program EERAD2 2 , which 
computes the four-jet production at NLO, we 
implemented the NNLO antenna subtraction 
method for the l/N 2 colour factor contribution 
to e + e~ — > 3j. EERAD2 already contains the five- 
parton and four-parton matrix elements relevant 
here, as well as the NLO-type subtraction terms. 

The implementation contains three channels, 
classified by their partonic multiplicity: (a) in 
the five-parton channel, we integrate da^ NLO — 
d- a NNLO'> 0°) m the four-parton channel, we inte- 
grate d<rJ^ LO — d<rj^£ ; (c) in the three-parton 

channel, we integrate dcrJ^ LO + daff NLO + 
V s 1 

d<T N pj L0 . The numerical integration over these 
channels is carried out by Monte Carlo methods. 

By construction, the integrands in the four- 
parton and three-parton channel are free of ex- 
plicit infrared poles. In the five-parton and four- 
parton channel, we tested the proper implementa- 
tion of the subtraction by generating trajectories 
of phase space points approaching a given single 
or double unresolved limit. Along these trajec- 
tories, we observe that the antenna subtraction 



terms converge locally towards the physical ma- 
trix elements, and that the cancellations among 
individual contributions to the subtraction terms 
take place as expected. Moreover, we checked 
the correctness of the subtraction by introducing 
a lower cut (slicing parameter) on the phase space 
variables, and observing that our results are in- 
dependent of this cut (provided it is chosen small 
enough). This behaviour indicates that the sub- 
traction terms ensure that the contribution of po- 
tentially singular regions of the final state phase 
space does not contribute to the numerical inte- 
grals, but is accounted for analytically. 

Finally, we noted in yQ that the infrared poles 
of the two-loop (including one-loop times one- 
loop) correction to 7* — > qqg are cancelled in 
all colour factors by a combination of integrated 
three-parton and four-parton antenna functions. 
This highly non-trivial cancellation clearly illus- 
trates that the antenna functions derived here 
correctly approximate QCD matrix elements in 
all infrared singular limits at NNLO. They also 
outline the structure of infrared cancellations in 
e + e~ — ► 3j at NNLO, and indicate the structure 
of the subtraction terms in all colour factors. 

5. Outlook 

In this talk, we presented a new method for the 
subtraction of infrared singularities in the calcu- 
lation of jet observables at NNLO. We introduced 
subtraction terms for double real radiation at tree 
level and single real radiation at one loop based on 
antenna functions. These antenna functions de- 
scribe the colour-ordered radiation of unresolved 
partons between a pair of hard (radiator) partons. 
All antenna functions at NLO and NNLO can be 
derived systematically from physical matrix el- 
ements. To demonstrate the application of our 
new method on a non-trivial example, we imple- 
mented the NNLO corrections to the subleading 
colour contribution to e + e~ — > 3 jets. 

An immediate application of the method pre- 
sented here is the calculation of the full NNLO 
corrections to e + e" — > 3 jets [201 - T ne antenna 
subtraction method can be further generalised to 
NNLO corrections to jet production in lepton- 
hadron or hadron-hadron collisions. In these 
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kincmatical situations, the subtraction terms are 
constructed using the same antenna functions, 
but in different phase space configurations: in- 
stead of the 1 — ► n decay kinematics consid- 
ered here, 2 — > n scattering kinematics are re- 
quired, which can also contain singular configu- 
rations due to single or double initial state ra- 
diation. These require new sets of integrated 
antenna functions, accounting for the different 
phase space configurations in these cases. 
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